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1 Introduction 

It has been recognized that the Birkhoff-Gustavson normalization |^] works effectively in study- 
ing various nonlinear dynamical systems. For example, when a two-degree-of-freedom Hamilto- 
nian system with a l:l-resonant equilibrium point is given, for example, the BG- normalization 
of its Hamiltonian around the equilibrium point provides an 'approximate' Hamiltonian system 
whose Hamiltonian is the truncation of the normalized Hamiltonian up to a finite degree: The 
approximate Hamiltonian system provides a good account of the surface of section with suffi- 
ciently small energies of the given system [^, |^. Such a good approximation thereby implies that 
to find the family of Hamiltonian systems sharing the same BG-normalization up to a finite- 
degree amounts to find a family of Hamiltonian systems admitting the surface of section similar 
to each other. The following question has been posed by the author as the inverse problem of 
the BG-normalization [Q, ^, ^, ^ : What kind of polynomial Hamiltonian can he brought into a 
given polynomial Hamiltonian in BG-normal form, ? 

The Birkhoff-Gustavson normalization to be dealt with this paper is outlined as follows (cf. 
1^, H]): Let R" x R" be the phase space endowed with the Cartesian coordinates (g,p) working 
as canonical coordinates. Let K{q,p) be a Hamiltonian function defined on R" x R" (or a 
certain domain of it) which admits the power-series expansion 

n oo 

Kiq,p) = E Yip] + 5i) + E Kkiq,p), (1) 

i=l fe=3 

around the origin of R" x R", where each Kk{q,p) (A; = 3,4, • • •) is a homogeneous polynomial 
of degree-A; in {q,p), and {i^j} non- vanishing constants. 

Remark 1 The convergent radius of the power series ^ may vanish ^; this happens to any 
K that is not analytic hut dijjerentiahle around the origin, for example. In such a case, the 
power series ^ is considered only in a formal sense. We will, however, often eliminate the 
word 'formal ' from such formal power series henceforth. 

The BG-normalization oi K{q,p) is made as follows: Let (^,??) be another canonical coordinates 
of R" X R" and let the power series (see Remark |l]) 

n oo 

Wiq,v)=J2'lj^^ + T.^kiq,v) (2) 

j=l fc=3 

be a second-type generating function M, a function in the old position variables q and the new 
momentum variables rj, associated with the canonical transformation 

dS dS 

{q,p) ^ {^,ri) with P=-Q- and C = -^, (3) 

where each Wk{q,ri) (k = 3,4, • • •) is the homogeneous polynomials of degree k in {q,rj). Note 
that the transformation (^) leaves the origin of R" x R" invariant. Through the transformation 
(^, the Hamiltonian K{q,p) is brought into a power series, say G{^,r]), subject to 

dW dW 

It is easy to see from (|l|), (|2|) and @ that G{S,,ri) takes the form 

n oo 

G{tv) = E T (^' + ^') +Y.Gk{tv), (5) 

j=l ^ fc=3 

where each Gk{£,, rj) (A; = 3, 4, • • •) is a homogeneous polynomial of degree-A in (,^, r]). 



Definition 1.1 The power series G{(,,r}) is said to be in the Birkhoff-Gustavson (BG-) normal 
form up to degree-r if G{^, r]) satisfies the Pois son- commuting relation, 

Gk{^,v)\ =0 ik = 3,---,r), (6) 




where {■,-}^,-q denotes the canonical Poisson bracket to the coordinates {S,,ri) (see JT^/j. 

The inverse problem of the BG-normahzation has been hence posed as follows |l], ^: For a 
given power series Hamiltonian G{^, r/) in the BG-normal form ^ with (^ r = ooj, identify 
all the possible power-series Hamiltonians whish share G{^, rj) as their BG -normalization. In 
contrast with the inverse prblem, we will refer to the problem of normalizing Hamiltonians into 
BG-normal form as the ordinary problem. 

Since elementary algebraic operations, differentiation, and integration of polynomials have 
to be repeated many times to solve the inverse problem, computer algebra is worth applying 
to the inverse problem; see [|, |, |n| for 'ANFER' and [|, g for 'GITA"^'. The procedure in 
ANFER is based deeply on the composition of canonical transformation, which is published in 
P] only for the two-degree-of- freedom in 1:1 resonance case. 

The aim of this paper is to present the procedure in ANFER in mathematical terminol- 
ogy. After the procedure, a new aspect of the Bertrand-Darboux integrability found through 
the inverse problem of certain perturbed oscillator Hamiltonians |Q is presented briefly as an 
application of ANFER. The contents of this paper is organized as follows. 

Section 2 sets up the inverse problem of the BG-normalization: For a better understanding of 
the inverse problem, the solution of the ordinary problem is given first. After that, the inverse 
problem is posed and solved. Section 3 is devoted to studying the composition of canonical 
transformations which provides a key of organizing ANFER. Those who have taken a course of 
analytical mechanics might think that the composition is well-known already. However, it seems 
that almost all the things known of are on the composition of infinitesimal transformations. 
Since the transformations dealt with here are not infinitesimal, section 3 is indeed important 
to organize ANFER. In section 4, the procedure in ANFER is described. Section 5 is for an 
application of ANFER: a new aspect of the Bertrand-Darboux integrability condition for certain 
the perturbed harmonic oscillators [Q] is presented. Section 6 is for the concluding remarks. 

2 The inverse problem of the BG-normalization 

The major part of this section is devoted to the inverse problem of the BG-normalization posed 
in H], p, m, ^. Before the inverse problem, however, we show the way solve the ordinary problem, 
which will promote a better understanding of the inverse problem. 

2.1 Solving the ordinary problem 

Let us start with equating the homogeneous-polynomial part of degree- /c (A; = 3, 4, • • •) in (^). 
Then equation (Q) is put into the series of equations, 

Gkiq,v) + iDg,^Wk)=Kk{q,r])+^kiq,r]) (A: = 3,4, • • •), (7) 

where Dg^^i is the differential operator, 

^ ^ f d d \ 



j=i 



dr]j dqj 



The $fc(g, ry) in (0) is the homogeneous polynomial of degree-A; in (g, r/) which is uniquely deter- 
mined by W3, ■ • • , Wk-i, K3, • • • , Kk-i, G3, • • • , Gk-i given: In particular, we have $3(9, 1]) = 
and 
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1 dW:^\ dK: 



'^4('^'^)=Eohe + 



j^i 1 2 V ^9j / 9pj 
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ai^s 1 (dW3\ dG3 



dqj 2 \ 977j / d^j 



(q^v) 



tl' <'> 



while <l>jt with larger A: would be of very complicated form to be described. 

To solve equation (0), the direct-sum decomposition induced by Dq^ri of the spaces of homo- 
geneous polynomials are of great use. Let us denote by Vk{q, if) the vector space of homogeneous 
polynomials of degree-A; in (g, rf) with real- valued coefficients [k = 0, 1, • • •). Since the differential 
operator Dq^ri acts linearly on each Vfc(g, 77), the action of D^^^ naturally induces the direct-sum 
decomposition, 

Vk{q, r?) = image Z?^ ker dJ^ (A: = 0, 1, ■ • •), (10) 

of Vfc(g,r/), where Dq^j) denotes the restriction, 

Df)i = D.Av.i,,,) (^ = 0'i'---)- (11) 

For Dq^^ and Dq.rj (A; = 3, 4, • • •), we have the following easy to be shown. 
Lemma 2.1 Equation ^) is equivalent to 

(^,,,(G,|(,,,)))(g,7?) = (Z)W(GfcI(,,,)))(<7,r?)=0 (A; = 3, • • • ,r). (12) 

Namely, 

Gfc|(,,^)GkerZ)(^^) (fc = 3,---,r). (13) 

Lemma 2.2 Let the coefficients, {i>j}j=i^...^n, 0/ X]?=i(^j/2)(p? + q]) in K{q,p) be said to be 
independent over Z (integers) if and only if 

n 

E ^J'^J = (^J e Z) ^ 7^. = (j = 1, • • • , n). (14) 

i=i 

holds true. Then, ker Dq^^ is spanned by the even degree polynomials, 

n 

Pm = ]^(^j + (IjY^^ ("^i • non-negative integer, j = 1, ■ • • ,n), (15) 

i=i 

if and only if {vj} are independent over Z (cf. W). In particular, if {vj} are independent over 
Z then kerD^^jj = {0} for every odd k. 

We proceed to solving (|^) now. According to (p^), let us decompose Kk{q,rj) and (^k{(l-,fl) 
{k = 3,4, •••) to be 

where 

KT'' (g, v) , ^T" {q, V) e image 2^ ^ , 

E:r(g,r?),c|>|-(g,^)GkerZ)W,. 



(k) I 1 (k) 

Since Gk{q,'r]) e kerDg,,j by Lemma |2.1| and since Dg^rjWk € image -Dg,,j , we obtain 

Gk{q,v)=Kt'{q,r]) + ^f'{q,v) (A; = 3,4, • • •), (18) 

as a solution of (0), where VFfc is chosen to be 

Wk{q, v) = (4S r' n(.) (^r" i(g,.) + ^r" )) (q^ v)- (19) 

\ ^' ' image _Llg_,, / 

What is crucial of (19) is that Wk{q, rj) € image Dg,,j (A: = 3, 4, • • •), which ensures the uniqueness 
of G(^, 7]): For a certain integer k> 3, let us consider the sum, W^ = VFK+(any polynomial in ker Dq^^), 
where G^ and W^ with k < k are given by ( p!^JT9| ). Even after such a modification, W,^ satisfy 
([T9|) still, which will lead another series of solutions of (0) with k > k. Therefore, under the 

restriction, W^ G image Dg,,j (A; = 3, 4, ■ ■ •), (1^) with (13) is said to be the unique solution of 
@. To summarize, the ordinary problem is denned as follows: 

Definition 2.3 (The ordinary problem) For a given Hamiltonian K{q,p) in power series 
(1), bring K{q,p) into the BG-normal form G{£,,r]) in power series (^ which satisfy ^^ and 
(6), where the second-type generating function W of the form (^ is chosen to satisfy ^ and 

t^fc(g,77)G image i5(^j (fe = 3,4, • • •)■ (20) 

Theorem 2.4 The BG-normal form G{^,rj) for the Hamiltonian K{q,p) is given by (^ with 
(11), where the second-type generating function W{q,r]) in power series (^ is chosen to satisfy 
(U) and (M). 



2.2 The inverse problem 

To pose the inverse problem appropriately, we start with looking the key equation, (^), of the 
ordinary problem into more detail from a viewpoint of canonical transformations. With W{q, rf), 
let us associate the inverse canonical transformation, 

(C,r?)^(g,p) with e = -^^^ and p=-^i^, (21) 

of (^), so that —W{q,vi) is regarded as a third-type generating function (Goldstein 1950), a 

function of the old momentum variables r] and the new position variables q. Equation (^) 

rewritten as 

d{-W) ^ d{-W) 

K{q,-^—)=Gi-^^,v), (22) 

is then combined with (p|) to show the following. 

Lemma 2.5 Let G{^,r]) of ^ be the BG-normal form for the Hamiltonian K{q,p) of ^, 
which satisfies (Op with a second-type generating function W{q,r]). The Hamiltonian K{q,p) is 
restored from G{^, r]) through the canonical transformation (pi) associated with the third-type 
generating function —W{q, rj). 

Now we can pose the inverse problem in the following way: Let the Hamiltonian H{q,p) be 
written in the form, 

1 n oo 

Hiq,p) = - E {P'j + l') + E Hkiq,p), (23) 

i=l fc=3 



where each Hk{q,p) (/c = 3,4, • • •) is a homogeneous polynomial of degree-k in {q,p). Further, 
let a third-type generating function S{q, rj) be written in the form, 

n oo 

Siq,v) = -J2^J^i~J2^kiq,v), (24) 

j=l k=3 

where each Sk{q,p) (fe = 3, 4, • • •) is a homogeneous polynomial of degree-fc in {q,p)- 

Definition 2.6 (The inverse problem) For a given BG-normal form, G{^,r]), in power se- 
ries ^), identify all the Hamiltonians H{q,p) in power series (^^ which satisfy 

Hiq,--) = Gi--,,), (25) 



where the third-type generating function S{q, r/) in power series ^4) is chosen to satisfy ^t 
and 

Sk{q,v) e image dJ^ (fc = 3,4, • • •)• (26) 

2.3 Solving the inverse problem 

We solve the inverse problem in the following way. On equating the homogeneous-polynomial 
part of degvee-k in (P5|), equation (|25|) is put into the series of equations. 



Hk{q,r,)-{Dg,^Sk){q,v)=Gk(.q,r])-^k{q,v) (A: = 3, 4, • • •), (27) 

where Dg^^ is given by (^). The ^k{q,v) is the homogeneous polynomial of degiee-k in (g, ??) 
determined uniquely by H3, ■ ■ •, H^-i, G3, ■ ■ •, Gk-i, S3, ■ ■ ■, Sk-i given. In particular, we have 
^3 ((7,7/) = and 



2 



1 dS3\ dH: 



^4(.,r,) = EUh^ + 



^■=1 y \dqj J dpj 



2 



9^3 I dS3\ dG. 



{q,rj) ^'i^ ^ V ^"^^ J ^^^ 




(28) 



while ^k with larger k would be of very complicated form to be described. Like in the ordinary 
problem, we solve (|27|) for Hk and Sk by using the direct-sum decomposition (|lO|) of Vfe(g,ry), 
the vector spaces of homogeneous polynomials of degree-A; (A; = 3, 4, • • •). Let us decompose H^ 
and ^k to be 

Hk{q,r^)=Hr''{q,^)+Hf'{q,^), 

^fe(g,7?)=rr"(g,r?) + M/^(g,r/), 

where 

Fr^°(<7,r/),v&r"(g,r?) G image Z^g, 

^r(g,^),^r(g,^)GkerZ?W, 
Then on equating ker Dg,,j-part in (p7[), H^'^^ is determined to be 



(29) 



(30) 



Hl<^^{q,r^)=Gk{q,r,)-^f^{q,r^). (31) 

Equating image L'g^^-part of ( p7| ) , we have 

H'T" [q, V) - iD,,,Sk) {q, V) = -^r" (Q, v)- (32) 



Since we have the pah' of unidentified polynomials, H^^'^'^" and Sk in (32), //™''*5° is not deter- 
mined uniquely in contrast with ffj^""" ; such a non- uniqueness is of the very nature of the inverse 
problem. Accordingly, equation (|3^) is solved to as 



H'^^''^'' {q,r]) E image D^^^ : chosen arbitrarily, 
-1 



(A: = 3, 4,. 



\ image Uq,, 

Thus we have the following. 



1(9.^) 



{q,v), 



(33) 
(34) 



Theorem 2.7 For a given BG-normal form G{^,r]) in power series ^, the solution H{q,p) 
of the inverse problem is given by (^^ subject to (31_) and p^J, where the third-type generating 
function S{q,r]) in (p3j is chosen to be (fl^ subject to (34)- 



Remark 2 We have another expression of and Sk and 7/™''^° equivalent to ( pl\, \33i , \34); 

H'r" (q, ri) = {Di'iSk) {q, v) + *r" («- v), (35) 

with 



Sk{q,rj) G image Z)^^ : chosen arbitrarily. 



(36) 



Since we usually pay much more interests in H than in S, in the inverse problem, it would 
better to take the expression (31, p^, ^2 



2.4 The degree-p ordinary and inverse problems 

From a practical point of view, we usually deal with the BG-normal form Hamiltonians not in 
power series but in polynomial. Indeed, as mentioned in section 1, when we utilize the BG- 
normalization to provide an approximate system for a given system, we truncate the normalized 
Hamiltonian up to a finite degree. Hence we naturally come to think of a 'finite-degree version' 
of both the ordinary and the inverse problems [|^, ^ . 

Definition 2.8 (The degree-p ordinary problem) For a given Hamiltonian K{q,p) of the 
form (Oj (possibly in polynomial form) and an integer p > 3, bring K(q,p) into the polynomial 
G{^, rj) of degree-p in BG-normal form which satisfy (Op up to degree-p, where the second-type 
generating function W{q, rj) in (0j is chosen to be the polynomial of degree-p subject to ^(\ ) with 
k = 3,---,p. 



Remark 3 In the papers /Q, ^, the degree, p, in Definition ^71 is restricted to be an even 
integers, because all the Vj 's are set to be 1 there, (cf. Lemmas 2A and 2A) However, since we 
are dealing with the general {vj} in this paper, the degree p has no restriction. 



Definition 2.9 (The degree-p inverse problem) For a given BG-normal form, G{^,rj), of 
degree-p with an integer p > 3, identify all the polynomial- Hamiltonians H{q,p) of degree-p 
which satisfy ^dj) up to degree-p, where the third-type generating function S{q, rj) is chosen to 
be the polynomial of degree-p subject to (p^ with k = 3,- ■ ■ , p. 



3 Composition of canonical transformations 

In Section 3, we have shown the way to solve the inverse problem which is theoretically complete. 
From a practical point of view, however, it is rather difficult to realize the solution in the present 
form, (31, 33, gj), since the calculations required will be highly combinatorial: Let us take the 1:1 
resonant two-degree-of-freedon case (n = 2, i^i = 1/2 = 2), for example. In order to calculate $fc 
(resp. ^fc), the polynomials G3, • • • , Gk-i, Sa, • • • , Sk-i, and Ks,---, Kk-i (resp. H^,---, Hk-i) 
have to be kept each of which is in V^ of combinatorially rising dimension, Yl,h=i ih) (/i-i)' where 
the symbol {') indicates the binomial coefficient. 

To get rid of such a difficulty, we break the needed calculations into a series of calculation 
of simpler form, which is realized in ANFER (see H, |2|, |ll|). This section is hence devoted to 
the mathematical basis for organizing ANFER. 



A key idea in writing-up ANFER is to realize (|3l|, ^, ^4D up to degiee-p by applying not 
the canonical transformation associated with S{q, rj) but the composition of canonical transfor- 
mations 

Th : {^^''-'\V^''-'^) ^ i&\v^^^) with (e(2),^(2)) = (^,^) (^r = 3,---,p), (37) 

associated with the third- type generating functions, 

5W(^W,^(^-i)) = _;f^^W^(.^-i)_5,(^W,^(/^-i)) (/, = 3,...,p), (38) 

i=i 

where Sh{C Kv ) is the homogeneous polynomial part of degree-Zi of the generating function 
S{q,ri) (see (H)) with (^W,r/('^-i)) in place of (5,7?). 

We start by showing the following on the generating function associated with the composition 
of a pair of canonical transformations. It should be pointed out that little is known explicitly 
of the generating function of the composition of non-infinitesimal canonical transformations like 
the following lemma, while well-known is that of infinitesimal ones (see Q, for example). 

Lemma 3.1 Let the function fw(uW ^v^^^^) (h = 1,2) be the third-type generating functions 
of canonical transformations, 

^;,:(n('^-i),T;(^-i))^(nW,t;('^)) (/i = l,2), (39) 

with 

^h-i) _ __df}^ ^h) _ _df^ (h-l2) (40) 

each of which leaves the origin o/R"xR" invariant. Then the composition, a2oai : {u''^',v^^') — > 
(u^'^' ,v^'^'), is generated by the third-type generating function, 

n 

/,,2(n(2),^(0)) = ^S,.i),. + /a)(S,.(0)) + /(2)(n(2),i)), (41) 

where u^^' and v^^' on the rhs of (0j are the C°° -functions of {u^'^' ,v^^') uniquely determined 
around {u^'^\v^^') = (0,0) to satisfy 



Proof: Since the canonical relation ( ^0|) yields the identities ||9[ 

n 

- Y.{uf'^^dvf-^^ + vfduf^) = ^(/^(uW,!;^^-!))) {h = 1,2), (43) 



we have 



n / " 

i=i Vi=i 



Further, the canonical relation (^0|) restricts {(n^'^-', w*-^^)}/i=o,i,2 to the inverse image, a "^(0, 0, 0, 0), 
of the C°°-map 

a:(nW,^W,n«,r;«,u(2),^(2))eR6 

l" (0) ^ 5/^'^ (1) ^ ^/^'' (1) ^ 5/^'^ (2) ^ 5/W \ (45) 

Applying the implicit function theorem jl^ to the map o", we obtain uniquely the functions, ■u, 
V, u and v, of {u^'^\v^^') subject to a{u,v^^\u,v,u^'^\v) = around (0,0). The differentiability 
comes from that of a due to the implicit function theorem. The substitution of (tt, v) into 
{u^^' ,v^^') in (44) thereby shows (|4l| ) with (p2|) . This ends the proof. 

We are now in a position to compo se t he series of canonical transformations {Th}h=i,--,p^ 



given by ( P7[) and ( |38D by using Lemma 3^. We show the following for the generating function 
of the composition, Tp o • • • o T3. 

Lemma 3.2 For any fixed integer p > 4, the composition, TpO ■ ■ ■ o t^, of the canonical trans- 
formations {Th}h=z,-,p given by 

C''~'' = -^i^''\v^'-'^ V^'' = -^(^''\v''-'') (/^ = 3,4,...). (46) 



with ( ^^ ) and (37) is associated with the third-type-3 generating function S^'''{^^''' ,7]^^^') defined 
recursively by 

5(3)(e(3),r?(2)) = 5(3)(^(3),,^(2)), (47) 

5W(^W,^(2)) 

n 

i=i 

where S^ ' are the third-type generating functions given by ^SSj), and ^("~^) and fj^'^^^> the C°°- 
functions of {£,^^' ,r]^'^') determined uniquely to satisfy 

around (^('*),r?(2)) = (0,0). 

Proof: We prove the lemma by induction: The starting case, (H^) with h = 4, immediately 



follows due to Lemma 3T. We move on to show the case of /i = r(> 5) in t urn under the 



assumption that (|4§| ) and (^) hold true for h = r — 1. Applying Lemma p.l| to the pair of 

9 



transformations, r^. with S^"^' and r^-i o • • • o T3 with 5^*" ^', we obtain S^^''{^^^' ,7]^'^') given by 
([4^ ) and (|49| ) with h = r as the third-type generating function of r,. o (Tr_i o • • • o rs). The 
differentiabihty of .^'''~^'* and fj^'^~^> follows from that of the generating functions S^'''~'^' and S^^' 
(see the proof of Lemma ^J). This ends the proof. 

Expanding S^ ' in terms of {C ,V )i '^^ have the following, a key theorem of ANFER. 

Theorem 3.3 The third-type generating functions, S^ '{^^ \r]^'^') (h = 4, 5, ■••j, admit the 
power series expansion, 



n h 



(h) J2), 



i=i 



A:=3 



(E-=i(ef +^f ))^/^ 



(50) 



0), 



where Sk{£,^^' ,r]^^') are the homogeneous part of degree-k in S{q,r]) (see ^^) with {^^^',r]^'^') in 
place of {q,rj). 

Proof: We show (|50|) by induction. Expanding ( ^9|) with /i = 4 in terms of (^^^^r/'^^), we have 



dSA 



|(3)=^(4) + i!:£l(^{4)^-(3)) 



5r/(3) 



^(3) = ^(2) + 



95; 



0^(3) 



Irie^'^^^^^)- 



(51) 
(52) 



Since the second term on the rhs of ( 51^) is of degree higher than two and since that of ( |5^ ) are 
of degree higher than one, ([5T|) and (32) are put together to provide a further expansion, 



i^'^=(^'^ + i^{&\v^'^) + o's{^-',V 



fj(3) ^ ^(2) ^ 



dS4_ 
9ry(3) 

dSs 



(4) „(2)l 



(e('\^(^)) + 03'(e('\r?(2^), 



(53) 



,(4)^ 



.(2)^ 



where 03 and 03 indicate the terms in (\/X]?=i(^i + ^J- )) of degree higher than three. 
Equations (p^J4^J53|) are then put together to yield 



5(4) (^(4) ^^(2)^ 

V I ;^(^) +— ^ 



(3) 



+ O3 



(5(4),^(2)) 



(2) ^gs 



ser 



+0;; 



E I ^f + ""^^ 



i=i 



H 



(3) 



+ 0Ur/f + 53(e(^) 



9^4 



(g(4),r,(2)) 






(4) 



V 



(2) , dSs 



da 



Or/(3) 



(5(4),^(2)) 






953 



(^(4) ,^(2)) 



9^(3) 



(^(4) ,^(2)) 



+ 0^ 



t(4)„(2) 
3=1 k=3 



E^fvf -EsU^^'^^v^'^) + o.i^^'\v^'^] 



(54) 
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which shows (^) with h = 4. We move on to show (|5y) with h = r > 5 in turn under the 
assumption that (pO|) with h = r — I holds true. In a similar way to get (|53|), we obtain 



e^"-^) = &^ + 7^(e('-\r?(^)) + o;_,(e«,ry(2)), 



9r?(^-i) 



={r-l) _ „(2) 



r-1 



a^. 



(55) 



^^- -=ri-+E^ii^'^^V^'^)+o'UiC^^\v 



k=3 



d^(k) 



(2)> 



where o[._i and o"_]^ indicate the terms in {\lY^^=i{Cj + ^i )) of degree higher than r — 1. 



Like (|5^, equations ( pq ) with h = r, ( ^8[ ) with /i = r — 1 and (|^) are thereby put together to 
yield 



5M(^(r)^^{2), 



E U, 



(r) 



as. 



j=i 



5?7! 



(r-l) 



r-1 



95, 



(^('■),,,(2)) 



-;.||.r^i:-,,., 



+ o'U 



{^('■),r;(2)) 



/ 



^/ M^ as, 



i=i 



V 



9?7 



(r-l) 



r-l 



+ o;_,hf) 



(^('■),,,(2)) 



+ E'5fc(^^'^ + 



95, 



fc=3 



9r/(^) 



(^('■),,,(2)) 



+ o',_i,r/j^^; 



i=i 



-•J 



+ o: 



r-l 



r-l 



+Sr{e\v^'^ +T. ""' 



fc=3 






+ o: 



{^{r),r,{2)) 



r-1^ 



E^r^f - E ^^(e^^^^^'^) + -^(e^''^^^'^)' 



+5r(^w,r?(^): 

j=l fc=3 

which proves our assertion. 

Corollary 3.4 S^^'{q,r]) coincides with S{q,ri) up to degree-h if expanded. 



(56) 



4 The procedure in ANFER 

In the previous section, we have studied the composition of canonical transformations necessary 
to organize ANFER. This section is devoted to the solving procedure in ANFER in which the 
composition of transformations is utilized effectively. 

4.1 The solution of the degree-p inverse problem 

We show that the solution of the degree-p inverse problem is obtained by applying the composi- 
tion Tp o ••• o r3 of canonical transformations studied in section. 3. Let us fix a BG-normal form 
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power series (possibly a polynomial) Hamiltonian G{^,r]) of the form (g). For G{^,rj) fixed, we 
define the power-series Hamiltonians H^ ' {^^ \r]^ ') (/i = 3, • • • ,p) by the reccurent formula 

H^^UTh = H^^-^^ (/i = 3,4,---) with F(2)(^,r/)=G(e,r/). (57) 

Namely, H^^' satisfy the equation 

(F('^)oT/,o...oT3)(e,r/)=F(2)(e,?7)(=G(e,r/)) (^ = 3,4,...), (58) 

which is put together with Lemma |3.2| to show 



^""K""-f!J^) = «'^'(-|S."'^') (" = 3.4.^0. (59) 



where S^'^\s!^^\ri^'^^) is given by (|8|) with (|9|) 



Let us recall the equation ( p7| ) equivalent to the defining equation ( |25| ) of the inverse problem. 
From (|27|), we see that only the homogeneous parts in 5 and H of degree lower than h, and those 
in G of degree lower than or equal to h (cf. (P8|)) are necessary to determine the homogeneous 
part Hh of degree- /i in H This observation is thereby put together with Corollary ^^ to show 
the following. 

Theorem 4.1 For a given BG-normal form power-series G{^,r]), the power-series Hamilto- 
nian H^ '{q,p) coincides up to degree-h (h = 3,4,-- ■) with the solution of the inverse problem, 
H{q,p). 

Corollary 4.2 The solution of the degree-p inverse problem is given by H^P'{q,p) truncated up 
to degree-p. 

4.2 The solving procedure in ANFER 

We are now in a position to present the procedure of solving equation (p7|) with h = 3,- ■ ■ p 
(equivalently, ( p7D up to degree-p) for the degree-p inverse problem. In ANFER, the series of 
equations ( |57D are solved as follows: To be more precise, what we solve is 

^"■'K"".-f|^) = '^"'-"(-^.,"'-") (ft = 3,...,p) (60) 

with 



equivalent to (|5^). For convenience, we will refer to the stage of solving (pOJ ) with h = r as the 
'Stage-r' hence force. Before starting the Stage-three, (pi]) is assumed to have been proceeded 
already. Further, it is convenient to express H^ '{^^ ' ,r]^ ') to be 

-. n „ „ oo 

^W(ew,^W) = i^(,f Vef VE^f (e^'^^^'^) (62) 

j=l fc=3 

like H{q,p) (see (23)), where each Hf, '(^^(^\r]^^>^ (^k = 3,4,- ■ •) is the homogeneous polynomial 
part of degree-k in {S,^'^' ,r]^^'). 
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[Stage-r] : 

At the Stage-r, ( pOI) with h = r is solved. Equating the homogeneous part of degree-k {k 

3, • • • , p) in (pOl) with h = r, we have the series of equation, 



^M(^M^^(.-i)) ^ ^(-i)(eM,^(-i)) (A: = 3, • . . ,r - 1), (63) 

^M(^M^^(r-i)) ^ (z)^,,,_^(,_,)5,) (e(^),77('-i)) +i/(^~i)(^M,ry(^-i)) {k = r), (64) 
^M(^M^^(r-i)) ^ ^(r-i)(^M^^(._i)) ^ 0M(^M^^(.-i)) (^ = ^ + 1, . . . ,^). (65) 

Here &^ l^^"^' ,r]^^^^') {k = r + 1, • • • ,p) in (|6^) are the homogeneous polynomials of degree-A; 
given by 



ei^)(eM,r/('-i)) 



E 

|a|=l 



a\ 



dSr 



drj 



(^(r)^^(r-l))^ 



d 



H, 



(r-l) 



9e 



Q^{r-l) ) ""fe-(r-2)|a| 



/ a A" „ 



(gM,,,(--i)) 



(66) 



{& 



)^^(._i))y vV9r?M; ^^fc-{r-2)iai 



{?('■) ,?y('-i)) 



where [(A; — 2)/(r — 2)] stands for the integer-part of {k — 2)/(r — 2), and a = (ai, • • • , a„) does 
the multi-index with nonnegative-integer valued components associated with the notations, 



a 



^I = E«i' 




a! = ai\ ■ ■ 


•"n!, 




j=i 










95,. \" /a5,.\"i 


/a5A"" 


fdSrY 


(dSrV^^ 


fdSr\ 


9ey UeJ 


V^Cn/ 


\dr] J 


\dili) 


\dVn/ 


(9\° /a \"i 


/ Q xan 


f dy 


f d \"i 


f d Y 


sd " UeJ 


U^nJ ' 


\dr]) 


Kdm) 


\dr]nJ 



(67) 



Remark 4 We do not have to deal with the homogeneous part in ( \6ul) of the degree higher than 



p taking the observation to (21) into account made above Theorem 4-1 



Since Theorem [4.1| and (|63| ) are put together to imply that each H^. with k < r has been 



(r) 



r(fc) 



identified already as Hj^ at the stage-fc, we do not have to deal with 
replacement of Hj^ ' to Hj^'^' {k < r) is made. 

We solve (M) and (pSi) in turn by using the decomposition. 



a lot: Only the 



kcr 



^«(e«,,(-i)) = ^(:)— (e«,,(-i))+^(:)-^(e«,,(-i) 



{k 



,p) 



ker 



0M(^M,^(-i)) = e('^)'-"^^°(CM,^(-i)) + er)'^^(CW,^(^-^)) {k = r + l,---,p) 



where 



^W--(^M,^(-i)),eW^-^^^^(CW,r?(-^)) G image Z5 



(k) 

^('■),,;('--l)! 



F 



«'^^(C«,,(^-^)),er)'^^^(^«,,(-^) G kerZ)(f) ^,(.-.- 



(68) 



(69) 
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Like in the way to solve (^) in subsection 2.3, (|64[) is solved to be 

^M'^-(^(r)^^(r-i)) ^ ^(-i)'^'='-(^W,^{-i)) (70) 

^^r)--s= (^(r) ^ j^(^-i))) e image D^^^, ^(,_i) : chosen arbitrarily, (71) 

„M l-f^r - -^r l(5M,,7('-i))J (4'ST ')• I'-^j 

We move on to ( |65| ) with k = s {> r): Since we have identified Hj[ s with k < s and S'r already 
before solving ( |65| ) with A: = s(> r), ©s can be identified completely by ( |66| ) with k = s. After 
the calculation of 0^ , the Hj^ with k = s{> r) is thereby identified by (|65|) . 

In ANFER, the solving procedure described above is realized in the symbolic computing 
language, REDUCES. 6 or later. Till to now, a prototype of ANFER for the 1:1 resonant two- 
degree-of- freedom systems (n = 2, z^i = 1^2 = 1) has b een written up by the author. The source- 




code of this prototype is available on the web page, http://yang. amp. i.kyoto-u. ac.jp/ 'uwano/ 







In closing this section, we wish to compare the procedures, (p7|)-(p4|) with (|6^)-(|7^). Through 
|27|)-(|3^ with k = r, S^, • • • , ^^-i, G3, ■ ■ ■ ,Gr, and H^, ■ ■ ■ ,Hr have to be kept. In contrast with 



(r—l) (r—1) (r) (r) 

this, Sr-i, H^ ,■■■ , Hp , and H^ , • • • , Hp have to be kept. Hence, the procedure given 
in this section would contributes to the memory-saving. 

5 Application 

Although only a prototype of ANFER exists for a 1:1 resonant case, ANFER has worked very 
effectively to find a new deep relation between the Bertrand-Darboux integrability condition 
(BDIC) (see [ |l3i |lq , |ig| ], for example) for the perturbed harmonic oscillators with homogeneous- 
cuhic potentials (PHOCPs) and for the oscillators with homogeneous- ^itariic potentials (PHO- 
QPs) |1[. The BDIC provides not only a sufficient condition for the integrability but also a 
separability, which will be reviewed briefly in Appendix. 

5.1 The degree-4 ordinary and inverse problems for PHOCP 

We start by solving the degree-4 ordinary problem for the PHOCP-Hamiltonian 

1 ^ 
^^^\q,p) = 1; Y.ip] + 9|) + Uiql + f2ql + fmql + Uql), (73) 

where (fh) are real-valued parameters. Note that if we choose (fh) to be 

/l = /3 = 0, /2 = 1, /4 = //, (74) 

J-^^'{q,p) becomes the well-known one-parameter Henon-Heiles Hamiltonian |^, |^. 

Through ANFER ([11 1), we obtain the BG-normal form Q of the following form as the 



solution of the degree-4 ordinary problem 



^(e,^) = ^(ciCi + C2C^ 



- ^(/i c?Ci + flckl) - ^(/i/3CiC2CiC2 + /2/4C1C2C1C2) 
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5 , , , ,27 7 , f f /- /-T'^ , f f A ^T"^ : f f a2- 



(/1/2C1 C1C2 + /1/2C1C2C1 + /3/4C1C2C; + /3/4C2C1C2) 
(/2/3Cl'ClC2 + /2/3C1C2C1 + /2/3C1C2C2 + /2/3CIC1C2) (75) 



24 

- g(/2ClC2ClC2 + /3C1C2C1C2) - ^(/2ClCl + /3C2C2) 

- lifkKl + fkUl + fkfcl + /Iclci) 

+ -^ififscKl + fifsclcl + f2hcKl + /2/4ClCi). 

We solve the degree-4 inverse problem for the BG-normal form Q in turn: By ANFER, we 
have the following polynomial of degree-4 as the solution pi ; 

1 ^ 
n{q,p) = -Y.{p]+q])+n3{q,p)+n^{q,p), (76) 

J=l 
with 

T~(-3iq,p) = a-izf + a2zfz2 + 032:1 ^2 + a4z| + a^zfzi + aQzfz2 

+ 07^122^1 + asZiZ2Z2 + agzl'Zi + aioZ2Z2 

+ aizf + 02^1^2 + 0321^2 + 04^ + «5-Zl^l + 06-^2^1 

+ 0721^1^2 + 08-^2^1^2 + agzizl + aiQZ2Z2, (77) 

and 

T^A^q^P) = Cizf + C2zfz2 + C3zlz2+CiZiZ^ 

+C5Z2 + CQzfzi + Crzfz2 + C8z'fz2Zi 

+Cc,zfz2'Z2 + Cio^l^i^i + CUZ1Z2Z2 + C12Z2Z1 + Ci3zf Z2 

+Cizf + C2zfz2 + C3zlzl + C4Z1Z2 

+C5z| + CQZizf + C7Z2ZI + C82:i^i^2 

+C9Z2^1^2 + Cl02:i^l^2 + CllZ2^1^2 + Cl2^;i^2 + Cl3^;2^ 

+8(a6062;i22Zl^2 + a9a9ZiZ2^iZ2 + 0506^1^2^1 

+00052:1^1^2 + 09aio2;|ziZ2 + 010092:1 2:2^2) 

/ 'y \ 

+6(0101^1^1 + 04042:2^2 + 0505^121 + 010010-^2^2) 
+4(01022:1^1^2 + 0309^1^2 + 0304^1^2^2 + 05032:12:2^122 

+O6O7Z1Z2 + 06082:iZ2^2 + 07092:1^1^2 + 070102:12:2^1^2) 

+4(0201^1^2^1 + 04O32IZ1Z2 + 0305^12:2^1^2 + aja^z^il 
+asaQZ2ZiZ2 + 0907^12:221 + 0190721^2^122 + 0903222^) 
8, _ __ _ __ , 

+ -(02022i222i22 + 030321222122) 

+2( — 06062i2i + 07072121 + O3O32222 — 09092222) 

+2(010321 22 + 02O42122 + 0507212122 — 05092122 

-0603212122 - aQO.iQz'l'z'l + 0703212221 

+07032^2122 - 0709212223 + 03010212222) 

, 2 2 2 2 2 2 2 

+2(03012221 + 04022221 + 0705212221 — 09052221 
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-asaeZiZ2zl - aioa^z^zl + asaTzlziZ2 
+asa7ZiZ2zl - aga7Z2ZiZ2 + a\i^a^zl^^2) 

+ o ("2a3-Z?^1^2 + a2«3-Zl-22^2 + «3«2-Zl^2^1 + a3a2z|ziZ2) 

+-{a2a2zjzl + aaasz^za) 

-T^Uhlzl + fhlA) - -l{flhziZ2ZlZ2 + f2f4ZlZ2ZiZ2) 

lb 4 

-ziflhzjziZ2 + flf2ZlZ2zl + hf4ZlZ2zl + /3/4 21^1^2) 
o 

0/2/3 / 2 2 2 2 \ 

-—{Z1Z1Z2 + Z1Z2Z1 + Z1Z2Z2 + Z2Z1Z2) 

--ifiziZ2ZiZ2 + fiziZ2ZiZ2) 

X.oQO 000 000 000\ 'J/OOO 000\ 

~o(/2^1^2 + /2^2^1 + /3^1^2 + /3^2^l) ~ T^(/2^l2:i + f?,Z2Z2) 

+^ifif3zM + hhzl-zl + f2Uz!4 + f2kzl-zl), (78) 

where a/i (/i = 1, • • • , 10) and q {I = 1, • • • , 13) are the complex-valued parameters chosen 
arbitrarily, and /^ (fc = 1, • • • ,4) the real- valued parameters in T^'^'{q) (see (|73|)). Namely, we 
have 46-degree-of freedom in the solution, 7^, of the inverse problem of the PHOCP if {fk) fixed. 
Note that if {ah), (q) are chosen to be 

oi = ^3 = as = as = ag = 0, 

1 2,^i (79) 

2a2 = 2a6 = 07 = -, 804 = aio = — , 

and 

Q = (£=!,•• -,13), (80) 

respectively, and (/fc) to be (74), TL becomes equal to the one-parameter Henon-Heiles Hamil- 



tonian. After (|7q)-(|7q), one might understand the necessity of computer algebra in the inverse 
problem. 

5.2 The BDIC for PHOCPs and PHOQPs 

We wish to find the condition for (a/i), (q) and {fj.) to bring Ti into the PHOQP-Hamiltonian 

1 ^ 
J'^^Xq.p) = 7J E(Pj + 9i) + (51^1 + 92qlq2 + QsqIqI + 94<?i<?i + 55^1), (81) 

where (g^) are real- values parameters. A straightforward calculation shows the following (see |j|] 
for detail) . 

Theorem 5.1 A PHOCP-Hamiltonian !F^^' shares its BG -normal form with a certain PHOQP- 
Hamiltonian J-^'^' up to degree-4: if and only if the PHOCP-Hamiltonian T^^' satisfies the BDIC, 

3(/i/3 + /2/4)-(/| + /|) = 0, (82) 
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for PHOCPs. Under ^^), the PHOQP-Hamiltonian ^'^^ sharing its BG-normal form with the 
PHOCP-Hamiltonian J-^"^' is equal to 

Q = ^E(pK'?|)-^(9/i+/l)'?t-y(3/i + /3)/2g?«2-^(/| + /|k?gi 
i=i 

If) ^ 

--^(3/4 + f2)f3qiql - ^(9/1 + fM, (83) 



where (fh) are subject to the BDIC ^^). 



We remark here that the BDIC (p2[) appears again in Appendix. 

We are now in a position to show the integrabiUty of the PHOQP with Q with (//i) subject 
to (^). It is easily seen that that J-^^' becomes Q under the substitution, 



5 10 5 

5i = -Y^(9/? + /l), 52 = -- ^(3/i + /3)/2, 53 = -3(/| + /|)> 

If) K 

54 = -y(3/4 + /2)/3, 55 = -Y^(9/| + /|). 



(84) 



A long but straightforward calculation shows that {gi) given by ( |84D with (|82D satisfy the BDIC 

9(72 + 4^3 - 24<7i(73 - 9^254 = 0, 

954 + 4^3 - 245355 - 9^254 = 0, (85) 

(52 + 54)53 - 6(5154 + 5255) = 0, 

for PHOQPs (see Appendix). To summarize, we have the following. 

Theorem 5.2 // a PHOCP and a PHOCP share the same BG-normal form up to degree-A, 
then both oscillators are integrable in the sense that they satisfy the BDIC. 

6 Concluding Remarks 

We have described the procedure in ANFER and the application of ANFER to the BDIC for 
the PHOCPs. In the following. Several remarks are made on ANFER itself and its application 

(1) (Improving the program): As pointed out at the end of section 4, the procedure in ANFER 
based on the composition of canonical transformation would save save the required-memory 
comparing with the procedure given in section 2.3. There would be much room for im- 
provement on this direction. It would be done in writing-up process, so that the author, 
not so familiar with writing programs, is trying to improve ANFER with several collab- 



orators |11|. Even if restricted to for PHOCPs, the inverse problem is expected to have 
rich related subjects listed below. Hence, those who are interested on improving ANFER 
will be welcome writing-up, because the improvement on ANFER is expected to advance 
various studies related with the inverse problem of the BG-normalization. 

(2) (the separability): By a simple calculation, we see that the PHOCP-Hamiltonian with 
/i = 73/3 = a + b, 72/3 = fi = a — b and the PHOQP-Hamiltonian with gi = 53/6 = 55 = 
— 5(a^ + 6^), 52 = 54 = 20(a^ + b"^) share the same BG-normal form up to degree-4, where 
o and b are real- valued parameters ffl. These oscillator Hamiltonians are well known to 
be separable in qi ±52 (see |1^), Theorem *** is hence understood to provide a significant 
relation between the separability of the PHOCP-Hamiltonians and that of PHOQP. Since 
T^'^' with (p3) and Q are thought to include several classes of Hamiltonians separable 
in several coordinate systems other than gi =b (72, the separability will be worth studying 
extensively from the BG-normalization viewpoint in future. 
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(3) ( qua ntum bifurcation): Since the perturbed oscillators referred to in theorems ^^ and 
|5.2| are integrable, their quantum spectra are expected to be obtained exactly. These 
oscillators are hence expected to provide good examples of the quantum bifurcation in the 
BG-normalized Hamiltonian systems [2^, 21, [2^, ^] to study whether or not the quantum 



bifurcation in the BG-normalized Hamiltonian system for those oscillators approximates 
the bifurcation in these oscillators in a good extent. 

(4) (integrability) : As is easily seen, the solution (|7^ , ^ , [78|) of the inverse problem for 
Q admits fifty real-valued parameters. We may hence expect to obtain other integrable 
systems, so-called the electromagnetic type [|^], for example. 

On closing this section, we wish to mention again of the role of computer algebra in the inverse 
problem: Without computer algebra, for example, in section 5, it would have been very difficult 



to find theorems 5.1 and 5.2, 
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7 Appendix: the BDIC 

In this Appendix, we review the BDIC briefly. We will restrict our attention here to the two- 
degree-of-freedom natural Hamiltonian systems associated with the Euclidean metric because 
the PHOCPs and the PHOQPs dealt with in section 5 are of such type. 

The theorem due to Bertrand and Darboux is stated as follows (see [^, |lj] , for example) . 

Theorem 7.1 (Bertrand-Darboux) Let F be a natural Hamiltonian of the form, 

I ^ 

i=i 

where V{q) a differentiable function in q. Then, the following three statements are equivalent 
for the Hamiltonian system, with F. 

(1) There exists a set of real-valued constants, {a, (3, (3' ,'y,'y') ^ (0,0,0,0,0), for which V{q) 
satisfies 

' d'^V d'^V\ 

{-2aqiq2 - (3 q2 - Pqi + t) 




^ o ("92 - 79? + /3g2 - /3'gi + 7') (^^) 

dV dV 

+ 75-(6ag2 + 3/3) - 7^(6agi + 3/3') = 0. 
dqi dq2 

(2) The Hamiltonian system with F admits an integral of m,otion quadratic in Tnomenta. 

(3) The Hamiltonian F is separable in either Cartesian, polar, parabolic or elliptic coordinates. 



Due to the statement (2) in theorem 7.1, a natural Hamiltonian system with F is always inte- 
grable if ( |87[ ) holds true. In this regard, we refer to (p7|) as the Bertrand-Darboux integrability 
condition (BDIC). 



For the PHOCPs and the PHOQPs, Yamaguchi and Nambu |14] have given a more explicit 
expression of the BDIC (q^) convenient for section 5: 
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Lemma 7.2 Let !F^^\q,p) (k = 3, 4J he the Hamiltonians of the form ( fz|j and ^8^) . 

(1) For the PHOCP with J^^^>(q,p), the BDIC ^8^ is equivalent to either of the following 
conditions, ^), ^8^, or (jg^; 

3(/i/3 + /2/4)-(/| + /|) = 0, (88) 

/i = 2/3, /2 = /4 = 0, (89) 

/4 = 2/2, /i = /a = 0. (90) 



(2) For the PHOQP with T^^'{q,p), the BDIC (SJ) is equivalent to either of the following 
conditions, ( |gl| j or ^9^); 

(g, = 2g,=29,, 
[92=94 = 0, 

9^2 + 4<?i - 24^153 - 9<7254 = 0, 

9gl + Agl - 24^395 - 95294 = 0, (92) 

, (52 + 54)53 - 6(5154 + 5255) = 0. 



The equations, (|8q ) and (|94 ), are referred to as the BDIC (^) for PHOCPs and ( |85|) for 
PHOQPs, respectively, in section 5. 



In closing Appendix, we wish to mention little more of the BDIC: As stated in Theorem 7.1, 
the BDIC (^) provide not only a necessary and sufficient condition for an existence of first 
integrals quadratic in momenta [15, ^] but also for the separability in either Cartesian, polar 



parabolic or elliptic coordinates |1^|. Indeed, the BDIC has been studied repeatedly from various 



viewpoints; the separation of variables [13, !(:[, the complete integrability ^\, so-called the direct 



method [|l^] and the renormalization of Hamiltonian equation [14[, for example (see also [19[ as 



an older reference and the references in the above-cited literature). 
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